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ABSTRACT 

Autocatalytic reaction fronts between reacted and unreacted species may propagate as solitary 
waves, namely at a constant front velocity and with a stationary concentration profile, resulting 
from a balance between molecular diffusion and chemical reaction. A velocity field in the support- 
ing medium may affect the propagation of such fronts through different phenomena: convection, 
diffusion enhancement, front shape changes... We report here on an experimental study and lat- 
tice BGK numerical simulations of the effect of an oscillating flow on the autocatalytic reaction 
between iodate and arsenous acid in a Hclc-Shaw cell. In the low frequency range covered by the 
experiments, the front behavior is controlled by the flow across the gap and can be reproduced 
with 2D numerical simulations. An analytical solution is also derived under the condition of a 
weak flow velocity, and is found to be in reasonable agreement with our data. 
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Introduction 

Interface motion and reaction front propagation occur in a number of different areas , including 
flame propagation in combustion population dynamics 00] and atmospheric chemistry (ozone 
hole) . An autocatalytic reaction front between two reacting species propagates as a solitary wave, 
that is, at a constant front velocity and with a stationary front profile [SJEj. These issues were 
addressed earlier on, but only a few cases are understood, such as the pioneering works of Fisher 
and Kohnogorov-Petrovskii-Piskunov 0] on a reaction-diffusion equation with a second-order 
kinetics ^ [7\ |S] . Although the effect of an underlying flow on a flame propagation has been 
extensively analyzed , the advective effect on the behavior of an autocatalytic front has been 
addressed only recently [HlEB^DEl- In this case, the evolution of the concentration of each 
chemical species is given by the Advection-Diffusion-Rcaction (ADR) equation: 

8(1 — > — 

— + U.~vC = D m AC + af(C) (1) 

where C is the normalized concentration of the (autocatalytic) reactant, U is the flow velocity, 
D m is the molecular diffusion coefficient and a is the reaction rate. 



In the absence of flow ( U = ) , the balance between diffusion and reaction leads to a solitary 
wave of constant velocity V x and width l x . For the autocatalytic Iodate-Arsenous Acid (IAA) 
reaction studied here, the kinetics is of the third order £Q, /(C) = C 2 (l — C), and the following 
ID solution of equation (JTJ is obtained pfl I13j: 

C "'"= i + .».U W . ■ '» = V^r- y *°V^r < 2 > 

where z is the direction of the front propagation. 

For a reaction propagating along the direction of a unidirectional stationary flow, U, two 
regimes have been described, depending on the ratio r\ = b/2l x , where b is the typical length scale 
transversely to the flow direction |14l 1151 IT?)] . In the eikonal regime, r\ 3> 1, the front propagates 
as a planar wave, at a velocity given by the sum of V x and of the algebraic maximum of the 
flow velocity (projected onto the direction of V x ), and takes the according stationary form. In 
the mixing regime, 77 <gc 1 , the interplay between diffusion and advection enhances the mixing of 
the chemical species and leads to an overall macroscopic diffusion known as Taylor hydrodynamic 
dispersion ^7]. As a result, the front moves faster. However, it is still described by equation 
in which the molecular diffusion coefficient D m has to be replaced by its effective macroscopic 
counterpart. 

The main idea of the present paper is to address the effect of an unsteady flow on the front 
propagation. We measure, experimentally and numerically, the velocity and width of a chemical 
front submitted to a time periodic flow, of period T = 1/ f = 2tt/lu. The question of the relevant 
time scale, to which the time scale of the flow, T, has to be compared, is discussed. 

We extend the theoretical work by Nolen and Xin et al. JJ5], who derived the time-averaged 
chemical front velocity in an oscillating flow, by analyzing the temporal evolution of the front 
velocity. We note that in the tracer case (without reaction), Chatwin ^H] and Smith JHj showed, 
using a Taylor- like approach |17j . that a pulsating flow results in an effective time dependent dif- 
fusion coefficient, the time-average of which is larger than the molecular diffusion coefficient |20| . 

In this paper, we study, experimentally and numerically, a third-order autocatalytic Iodate- 
Arsenous Acid (IAA) reaction submitted to a pulsative flow. In section 1, we present the experi- 
mental set-up and the measurements obtained using a large set of frequencies and amplitudes of 
oscillations. In section 2, we compare the experimental results with 2D numerical simulations and 
we investigate a wider range of parameters with additional simulations. In the last section, we 
extend the theoretical result by Nolen and Xin ^2 to derive the temporal variations of the front 
velocity. 
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Experimental set-up and data 

We use the third-order autocatalytic Iodate-Arsenous Acid (IAA) reaction. The reaction front is 
detected with starch, leading to a dark blue signature of the transient iodine as the reaction occurs 
In the absence of flow, a reaction front travels like a solitary wave, with a constant 
velocity V£ xp ~ 20 fim/s and with a stationary concentration profile of width l exp ~ 100 /im. 
We study the front propagation in a Hele-Shaw (HS) cell of cross-section b x h = 0.4 x 8 mm 2 
(along x and y directions, respectively). The unidirectional (along z direction) oscillating flow is 
imposed at the bottom of a vertical HS cell, from a reservoir filled with unreacted species. This 
revervoir is closed with a thin elastic membrane, pressed in its middle by a rigid rod fixed at the 
center of a loudspeaker. Consequently, a displacement of a given volume of liquid in the reservoir 
induces a displacement of the same volume of liquid in the HS cell. The y — z plane of the HS 
cell is enlightened from behind and recorded with a CCD camera. The amplitude A and the 
pulsation lo = 2nf of the oscillating flow are imposed by the controlled sine tension applied to 
the loudspeaker, and measured in situ, from the recorded displacement of the air/liquid interface 
at the top of the partially filled HS cell. This displacement follows the expected A sin (uit) time 
dependence. Due to the constraint of our device, the imposed amplitude and frequency of the 
flow displacement are in the ranges A <E [0.07 mm, 1.7 mm) and / <G [0.01 Hz, 0.2 Hz] and the 
maximum velocity of the flow in the cell is roughly Um = Auj. The oscillating flow field in the HS 
cell is of the form (see Appendix for the full expression): 

U(x, y, t) = U M $te[f(x, y)} sin (wt) (3) 

The shape of the velocity profile depends drastically on the viscous penetration length l v = y/is/oj 
|21j . If l v is large compared to the cell thickness b (low frequency), the flow variations are slow 
enough for the steady state to be established. The resulting "oscillating stationary" velocity profile 
is parabolic in the gap and flat along the width h of the cell except in the vicinity of the side walls, 
in a layer of thickness b Conversely, for l v « b (high frequency), the fluid has not enough 
time to feel the effects of the solid boundaries and the velocity profile is flat over the whole cross- 
section b x h, except in the vicinity of each wall, in a layer of thickness l v . Figure ^ is a sketch of 
such an effect. For our dilute aqueous solutions of viscosity v w 10 -6 m 2 .s~ 1 and in our frequency 
range, the penetration length l u , which ranges between 0.8 mm and 4 mm, is larger than the cell 
thickness b = 0.4 mm. 

Hence, in most of our experiments, the "stationary velocity profile" is instantaneously reached, 
parabolic Poiseuille-like across the gap b, and almost invariant along the y direction (except in a 
layer of thickness l v close to the boundaries) . 

Figure El displays snapshots of a typical experiment: We do observe a front slightly deformed, 
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Figure 1: Sketch of the velocity profile in the gap of a Hele-Shaw cell for low frequency (LF, left) and 
high frequency (HF, right). 

propagating up and down (oscillating), with a downward averaged displacement, from the burnt 
product of the reaction to the fresh reactant. 

From this movie, the front is tracked and its location is plotted as a function of time. The 
so-obtained figure clearly shows the oscillation of the front position at roughly the imposed 
frequency and an overall drift of the front. 

The measurement of this drift in time leads to the time-averaged front velocity (V^ xp ). Figure 
0]displays (Vf xp ), normalized by V£ xp , versus the amplitude of the time-periodic flow field U = 
2Um/3 = 2Auj/3 (where 2/3 is the ratio of the gap-averaged velocity to the maximum one of a 2D 
gap Poiscuillc profile), also normalized by V£ xp . The increase of (Vj Xp ) with U is almost linear, 
with a slope slightly larger than 0.1. This demonstrates that the propagation velocity of a reaction 
front can be enhanced by a null in average, laminar flow. Moreover, as the mean advection in 
this time-periodic flow is zero, this effect comes clearly from some non-linear interplay. It could 
be attributed to the enhancement of the mixing due to the presence of the flow. 

As mentioned above, it is seen from the instantaneous velocity curve (Vj Xp (t), figure that 
the front velocity oscillates at the frequency of the flow. However, due to the experimental noise, 
it is difficult to obtain further information from this curve. 

We also noticed from our observation of the experimental movies that the width of the colored 
front, L(t), was likely to oscillate at a frequency twice that of the excitation, which, unfortunately, 
is not obvious on the static pictures (figure 2). We note that this feature could support the 
description of the front thickness in the framework of an effective diffusion of coefficient D, as the 
latter is expected to be insensitive to the flow direction and to depend only on the flow intensity 
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Figure 2: Time evolution of a chemical front in a pulsative flow field of amplitude A — 0.55 mm and 
period T = 50 s. Time increases from left to right and from top to bottom. Two images are separated by 
T/4 time intervals. The distance between two dashes is 1 mm. 

as D cx U 2 m3 (which here oscillates at 2/). 

To test this empirical observation, we measured the width L(t) of the dark blue ribbon. As 
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this ribbon corresponds to the presence of the transient iodine, L(t) is a qualitative measure of 
the chemical front width, but gives however the right time behavior. A classical Fourier analysis 
of L(t) was tried, but, due to the large amount of noise, it did not provide any reliable frequency 
dependence. Therefore, we used the more sensitive micro-Doppler method (see |2.'{| 121] and the 
references therein) which analyzes an instantaneous signal frequency The so-obtained oscillation 
frequencies /' of the width L(t) versus the imposed ones /, are displayed in figurc[SJ They collapse 
onto the straight line /' = 2/, which supports the contention that the front width oscillates at 
twice the frequency of the flow. 

In order to have further insight into the instantaneous features of the propagating front, we 
used numerical simulations, which are less noisy and give access to the behavior in the gap of the 
cell. 
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Figure 4: Normalized drift velocity of the chemical front (Vf Xp ) /V£ XF versus the normalized flow intensity 

U/V* xp . 
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Figure 5: Experimental oscillation frequency /' (Hz) of the front width versus the imposed oscillation 
frequency / of the flow (Hz). •: Experimental values obtained with the micro-Doppler algorithm, — : 
/' = 2/. 

2-D Numerical simulations 



Assuming a third-order autocatalytic reaction kinetics for the IAA reaction ^ |SJ I13| and a unidi- 
rectional flow U(x,y,t) in the z direction, equation Q reads: 
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dC TT dC _ fd 2 C d 2 C\ 



dt dz \ dx 2 dz 2 

where C is the concentration of the (autocatalytic) rcactant iodide, normalized by the initial con- 
centration of iodate, D m is the molecular diffusion coefficient and a is the kinetic rate coefficient 
of the reaction. The solution of equation (0J is obtained using a lattice Bhatnagar-Gross-Krook 
(BGK) method, shown to be efficient in similar contexts ^3 123 EH ■ The full 3D periodic 
velocity field U(x,y,t), in a Hele-Shaw cell, has been derived analytically in the Appendix (equa- 
tion H26|) ~). As mentioned above, the oscillating flow field does not vary much along the y direction 
(except in a boundary layer of the order of the gap thickness b) , and the velocity field, away from 
the side walls, is basically a 2D flow field, U(x, t), given by: 



U(x, t) = u M m,e 



^ cos (kx) 
~ cos(fc6/2) 



iuj( 



(5) 



where k = J ^ is a complex wave number. Note that for small frequency, lo -C vjb 2 (i.e. 
b l v ), equation (JSJ) reduces to an oscillating Poiseuille flow: U(x,t) « Um(1 — Ax 2 /b 2 ) sin (cut). 
The analytic flow field (equation (JSJ) is used in equation for the 2D simulation of the ADR 
equation by the lattice BGK method. 

In order to compare the results of the numerical simulations with the experiments, we used 
the same non dimensional quantities, namely b/l x xp (= 4), Um/V£ xp and the Schmidt number 
Sc = v/D m (= 500) which compares the viscous and mass diffusivities. The simulations were 
performed on a lattice of length A^, ranging between 2000 and 6000 nodes, and of constant width 
N x = 40 nodes during 2 x 10 5 to 4 x I0 6 iterations. The above experimental value of b/l^ p gives 
a numerical chemical length l x = 10 = y ^^- We chose D m = 5.10~ 3 and a = 10~ 4 , which 
sets the front velocity V x = \j aD ^ m = 5.10" 4 and the kinematic viscosity v — D m Sc = 2.5. The 
varying parameters in the simulations are the amplitude A and the frequency / of the imposed 
oscillating flow field. A typical movie of a numerical simulation is displayed in figure El 

It is seen from these movies that the front oscillates and travels from the burnt product to 
the fresh reactant. The mean concentration profiles are obtained by averaging along the lattice 
width, and analyzed along the same line as in the experiments. Figure0shows the time evolution 
of the displacement of the iso-concentration C = 0.5, obtained in the 2D simulations and in 
the experiments: The agreement between the two supports the contention that in our frequency 
range, the dynamics of the front is governed only by the variations of the velocity field in the 
gap (b = 0.4 mm), and that the (large) transverse extent of the plane of the experimental cells 
(h = 8 mm) plays no role. 

Figure [S] displays the resulting normalized drift front velocity (Vf) /V x versus the normalized 
flow velocity U /V x . The two sets of data, obtained experimentally and numerically, are also in 
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Figure 6: Numerical simulation of the front displacement (obtained with A — 10 and / = 2.5 x 10~ 6 ). 
The product of the reaction is in dark and the reactant is in grey. From top to bottom, time increases by 
200000 time steps (= 1/2 period). The lattice dimensions are 40 x 4000 (note that the aspect ratio of the 
pictures is not 1). 

good agreement. This leads us to analyze in more details the dynamics of the front in the gap, 
with the help of the 2D numerical simulations. 

The theoretical front width given in can be obtained from the second-order moment of the 
derivative of the concentration profile 16 : 



The use of the discrete version of © allows us to estimate the front width If. Figure displays 
the time evolution of If. After a transient time of the order of one flow period, the front width 
I f (t) becomes periodic and oscillates with a frequency twice that of the imposed flow, in agreement 
with the experimental result given in figure [S] 

For a stationary laminar flow, we showed in |lrj| that, under mixing regime conditions, the 
velocity and width of the reaction front result from a Taylor-like diffusion process, with an effective 
diffusion coefficient D e ff such that Vf = \/ aD e ff/2 and If = y/2D e f / / 'a, leading to Vf /If = a/2. 
An easy way to test the relevance of this effective diffusion description to the present case is to 
measure the ratio of the time-averaged values (Vf) / '{If). Figure ITU1 clearly shows that the relation 
between (Vf) and (If) is not linear. This allows us to discard the above description. 

The information given by both the experiments and the numerical simulations may be sum- 
marized as follows: In the presence of a periodic flow, the propagation of a chemical front in a 
Hele-Shaw cell is governed by the velocity profile in the gap. The front position drifts in the natu- 
ral direction of the chemical wave and undergoes oscillations at the frequency of the flow whereas 
the front width oscillates at twice this frequency. Finally, the front behavior cannot be described 
in terms of a Taylor-like effective diffusion. 

In the next section, we derive a model, under the assumption of weak flow velocity, and discuss 
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Figure 7: Time evolution of the front position. — : Experimental data, : Numerical simulations. The 

experimental values (A, f) (A in mm and / in Hz) are: (1.12, 0.01), (0.55, 0.01), (0.14, 0.08), (0.27, 0.04), 
(0.1,0.06), (0.14,0.04) from top left and clockwise. 

the above findings in the light of our theoretical approach. 



Theoretical determination of the front velocity 

We have derived in a previous work |16j . the velocity of a stationary reaction front, using a 
small parameter expansion method. The two parameters were the reduced gap thickness b/l x (or 
i] = b/2l x ) and the normalized flow velocity e = U /V x , with U = 2U~m/3 in 2D. We shall use 
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Figure 8: Normalized drift velocity (Vf)/V x versus the normalized flow velocity U/V x . •: Experimental 
data, □: Numerical simulations. 
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Figure 9: Time evolution of the normalized effective front width lf/l x for the numerical simulation A — 1 
and / = 10 -5 . The period T of the oscillating flow is indicated for comparison. 



the same method in order to derive the instantaneous front velocity and width. We note that this 
extends the work by Nolen and Xin J2| on the drift front velocity. 
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Figure 10: Normalized front velocity {Vf)/V x versus normalized front width (If) /l x - The line through the 
data (best quadratic fit) indicates that {Vf)/V x is not proportional to (lf)/l x - Dashed line : {Vf)/V x = 
(lf)/lx- 



The ADR equation in a 2D, unidirectional flow along z writes: 

dC rr dC o f d ' 2 ° d2 °\ rtfi rto 



(7) 



The 2D flow is imposed in a gap of size extension b. We assume a frequency small enough 
(l v = yjvjuj 3> b) to have an oscillating Poiseuille flow: 



4x 2 

U(x,t) = U u [l--p-)sm(wt) 



We also assume that the flow velocity is small compared to V x , 

e = \U M \/V X « 1 



(8) 



(9) 



and that, accordingly, the change in the front velocity due to the presence of the flow is small 
compared to V x , (F/(t) - V x ) < V x . 

Moreover, we assume that the concentration field is nearly uniform along the transverse x 
direction. Note that, for a passive tracer, this hypothesis is fulfilled when the Peclet number 
Pe = \U\b/D m is smaller than L/b, where L is the typical advection length (condition for Taylor- 
diffusion regime). In the presence of a reaction, this condition becomes erj 2 <C 1 EHj- We recall 
first the approach of Nolen and Xin |12| . in our notations. In the moving frame (s = z — Vft), 
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under the assumption © of a weak flow velocity (e <C 1), the concentration C(s,x,t) and the 
velocities U and Vf can be expanded in powers of e, as follows: 



C = C {s)+eC 1 (s,x,t) (10) 
U = eUx(x,t) (11) 
V f = V x +eV 1 (t)+e 2 V 2 (t) + ... (12) 

where C (s) = 1/(1 + e s /^) is the mean concentration profile (averaged over the gap and the time) 
and Ci(s, x, t) denotes deviations from the mean. Using the space and time Fourier decomposition 
of the flow velocity field, U\ = ^2 n b n e tknX+tuJt , for a monochromatic (to) velocity field (where 
k n = 2m: /b is the decomposition wave vector for a gap of width 6), and expanding the ADR 
equation in the moving frame, Nolen and Xin derived the time-averaged, drift velocity: 

(Vf) =V X (1+ e V2), with 7 = ^ £ (13) 

where fi = 47r b f '" is a characteristic frequency, proportional to the inverse of the typical diffusion 
time across the gap b. For the Poiseuille flow used here, one finds b n = (— l)™3/7r 2 n 2 . 

In this paper, we are interested in the time dependence of the front velocity, in the presence of 
a sine flow velocity (a single temporal mode, u>, in the Fourier decomposition of Ui). After some 
calculations and with the necessary assumption <C d d '^ 1 1 , we find (see [57] for details): 

=f(t) = l+£cos(wi)+eV[7i(l + cos(2c^)) + 72sin(2^)] (14) 

with 

71 = ^S» 2 (15) 

72 = ^S^nra (16) 

The equations Q14|). (|15|l and (|16|1 display no term involving u>/a, which demonstrates that the 
reaction kinetics does not play any role at this order. At first order in e the front velocity is the 
algebraic sum of the gap-averaged flow velocity and of the chemical wave velocity in the absence of 
flow. Note that this result is surprisingly similar to the one obtained for the front displacement in 
constant flows in the mixing regime [141 116| . The leading order pulsative flow contribution causes 
the front to oscillate with the flow frequency. At second order in e, the contributions to the front 
velocity are a constant one and one at twice the frequency of the flow field. Note that for u> = 0, 
we recover previous results: p£ w 1 + e + yj^e 2 ?? 2 1141 Averaging expression l|14|) over one 
period leads to the drift front velocity (fH5|l (Vj heo }/V x = 1 + E 2 rj 2r )i. We note that the amplitude 
K v of the velocity oscillation is mainly dominated by the order one term: K v ~ e. Figure 1111 
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shows that the theoretical results and the numerical simulation data for the drift velocity and its 
amplitude of oscillation are in reasonable agreement for small e (e < 1). 
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Figure 11: Left: Normalized drift front velocity {Vf)/V x versus e = \Um\/V x . □: Numerical simulation 
data, X: Theoretical predictions in the limit of small e. 

Right: Amplitude K v of the oscillation of Vf(t)/V x versus e. □: Simulations results, X: Theoretical 
results. 

The insets give the behavior of (V/)/V x and K v over a wider range of e. 



Conclusion 

In this paper, we have analyzed the influence of a time-periodic flow on the behavior of an autocat- 
alytic reaction front. The numerical simulations arc in reasonable agreement with the experiments 
and show that the front dynamics is controlled by the flow velocity in the gap of a Hclc-Shaw cell 
in the range of the parameters explored. We have shown that, contrary to the case of a stationary 
laminar flow, a Taylor-like approach cannot account for the (timc-avcraged) width and velocity 
of the front. The instantaneous front velocity has been derived theoretically, in the weak flow 
velocity regime. It is found to be in reasonable agreement with the numerical simulation results 
for \U M \/V X < 1. 
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Appendix : Determination of the oscillating velocity profile 
in a Hele-Shaw cell 



We consider a laminar flow, unidirectional along the z direction, in a Hele-Shaw cell of cross- 
section b x h (b <C h) in the x and y directions, respectively. The Navier-Stokes equation in the z 
direction writes: 



dU 
dt 



d 2 U d 2 U\ 1 dP 



" (17) 

9a; 2 dy 2 J p dz 

where a harmonic pressure gradient oscillating at frequency / = lu/2tt is imposed, and the flow 
velocity U(x,y,t) satisfies the boundary conditions: 



U(x,±^,t) = U(±^,y,t) = 



(18) 



In order to derive the analytical expression of U, we extend the method used by Gondret et 
al |221 to calculate the stationary flow field. First, we consider the flow field, written U* (x)e~ luJt , 
between two infinite planes (at x = ±| as in our simulations), which satisfies: 



~2 

d 2 u* 



l ^U* - — — 
dx 2 v pv dz 



leading to |21| : 



vjjp oz 



cos (kx) 
cos(f) 



(19) 



(20) 



where k = J ^ = (1 + i)-v/^- Following [22], we write now the solution of the full problem as: 



U(x,y,t) = {U*(x) + U**(x,y)}> 



where U**(x,y) must now satisfy: 



d 2 U** d 2 U** 



dx 2 



dy 2 



— U = 



with the boundary conditions: 



U**(±-,y) 
U**(x,±~) 



-U*(x) 



The solution of (|22|1 is found assuming a Fourier decomposition of the form: 

oc 

U**(x, y) = ^A n cos (k xn x) cos (k yn y) 



(21) 

(22) 

(23) 
(24) 

(25) 



where k xn = (2n — l)iv/b, in order to satisfy the boundary condition 1(24(1 . After some calculations, 
we obtain the solution of 1)17(1: 



U(x,y,t) 



dP 



up dz 



cos (kx) 



cos(M) 



n fc 2 cos (k xn x) cos (fcynj/) 
(2n - l)7rfc 2 n cos (fc^f) 



(26) 
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where ky n = ^ — k^ n . The velocity U(x,y,t) induced by an oscillating pressure gradient has a 
phase shift and a modulus which depend on space. We note that the maximum modulus Um occurs 
in the middle of the cell (x — y — 0). In Expression J3|JJ), the first two terms correspond to a 2D 
oscillating flow between two parallel boundaries at distance b (U*) [51], and the third one accounts 
for the finite size of the cell in the larger direction h ([/**). We note that at high frequency, a 
plug flow takes place in the section of the cell, except in a thin viscous layer l v = \fvjui close to 
the boundaries (y = ±h/2 and x = ±6/2). On the opposite, in the low frequency regime [l v ^> b), 
the oscillating flow has the same shape as a static one given in |22| . Hence, except in a thin layer 
(of thickness ~ b <C h) close to the side boundaries (y = ±/i/2), the flow is a parabolic Poiseuille 
flow of the form: 



References 

[1] S. K. Scott, Oxford University Press, Oxford (GB) (1994). 

[2] P. Clavin, Prog. Energy. Combust. Sci. 11, 1 (1985). 

[3] R. A. Fisher, Pro. Annu. Symp. Eugenics. Soc. 7, 355 (1937). 

[4] A. N. Kolmogorov, I. G. Petrovskii and N. S. Piskunov, Moscow Univ. Math. Bull. (Engl. 
Transl.) 1, 1 (1937). 

[5] P. A. Epik and N. S. Shub, Dokl. Akad. Nauk SSSR 100, 503 (1955). 

[6] A. Hanna, A. Saul and K. Showalter, J. Am. Chem. Soc 104, 3838 (1982). 

[7] Ya. B. Zeldovitch and D. A. Franck-Kamenetskii, Actu. Phys. USSR. 9, 341 (1938). 

[8] U. Ebert and W. van Saarloos, Physica D 146, 1 (2000). 

[9] G. Papanicolaou and J. Xin, J. Stat. Phys. 63, 915 (1991). 

[10] B. Audoly, H. Berestycki and Y. Pomeau, C. R. Acad. Sc. Pans, Series lib 328, 255 (2000). 

[11] M. Abel, A. Celani, D. Vergni and A. Vulpiani, Phys. Rev. E. 64, 046307 (2001). 

[12] J. Nolcn and J. Xin, SIAM J. Multiscale Modeling and Simulation 1, 554 (2003). 

[13] M. Bockmann and S. C. Miillcr, Phys. Rev. Lett. 85, 2506 (2000). 

[14] B. F. Edwards, Phys. Rev. Lett. 89, 104501 (2002). 

[15] M. Leconte J. Martin, N. Rakotomalala and D. Salin, Phys. Rev Lett 90, 128302 (2003). 




(27) 



16 



[16] M. Leconte J. Martin. N. Rakotomalala, D. Salin and Y. C. Yortsos, J. Chem. Phys. 120, 
7314 (2004). 

[17] G.I. Taylor, Proc. Roy. Soc. Lond. A 219, 186 (1953). G. I. Taylor, Proc. Roy. Soc. B 67, 
857 (1954). 

[18] P. C. Chatwin, J. Fluid. Mech. 71, 513 (1975). 
[19] R. Smith, J. Fluid. Mech. 114, 379 (1982). 
[20] E. J. Watson, J. Fluid. Mech. 133, 233 (1983). 

[21] L. Landau and E. Lifchitz, "Physique thcorique, mecanique des fluides" Mir edition, (1989). 

[22] P. Gondret, N. Rakotomalala, M. Rabaud, D. Salin and P. Watzky, Phys. Fluids. 9, 1841 
(1997). 

[23] I. Rcnhorn, C. Karlsson, D. Lctalick, M. Millnert and R. Rutgers, SPIE 2472, 23 (1995). 
[24] M. Lavicllc and C. Lcvy-Leduc, to appear in IEEE Transaction On Signal Processing, (2005). 
[25] J. Martin, N. Rakotomalala, D. Salin and M. Bockmann, Phys. Rev E, 65 051605 (2002). 
[26] N. Rakotomalala, D. Salin and P. Watzky, J. Fluid. Mech. 338, 277 (1997). 
[27] M. Leconte, "Propagation de front de reaction sous ecoulement", These de doctorat, (2004). 



17 



